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Abstract 

Spherical clusters of SU{2) BPS monopoles are investigated here. 
A large class of monopole solutions is found using an abelian approx- 
imation, where the clusters are spherically symmetric, although ex- 
act solutions cannot have this symmetry precisely. Monopole clusters 
generalise the Bolognesi magnetic bag solution of the same charge, 
but they are always larger. Selected density profiles give structures 
analogous to planets of uniform density, and galaxies with a density 
decaying as the inverse square of the distance from the centre. The 
Bolognesi bag itself has features analogous to a black hole, and this 
analogy between monopole clusters and astrophysical objects with or 
without black holes in their central region is developed further. It is 
also shown that certain exact, platonic monopoles of small charge have 
sizes and other features consistent with what is expected for magnetic 
bags. 
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1 Introduction 



SU(2) Yang-Mills gauge theory with an adjoint scalar Higgs field is a classic 
example of a nonlinear field theory in three space dimensions with smooth, 
topological soliton solutions. The solitons are magnetic monopoles [HEIE]- 
SU(2) magnetic monopoles were first discovered by 't Hooft and Polyakov, 
but the special case pinpointed by Bogomolny, Prasad and Sommerfield 
(BPS) has received particular attention from physicists and mathematicians. 
Static BPS monopoles obey a first order partial differential equation known 
as the Bogomolny equation, B = D$, where B is the Yang-Mills magnetic 
field and D$ is the covariant gradient of the Higgs field $. In addition, there 
is the boundary condition that at spatial infinity, $ is of unit magnitude; in 
other words, $ has unit vacuum expectation value. This Bogomolny equation 
is a three-dimensional analogue of the four- dimensional self-dual Yang-Mills 
equation, whose solutions are instantons. It is not easy to solve explicitly. 

Monopole solutions are characterised by a topological charge A^, a pos- 
itive integer that can be identified with the number of basic monopoles, so 
solutions with charge are A^-monopole configurations. The basic monopole 
of unit charge is spherically symmetric and has core size 1. Outside the core, 
the fields abelianize and resemble those of a Dirac monopole, with a magnetic 
Coulomb tail, and in addition a long-range Higgs field tail. 

A^-monopole solutions of the Bogomolny equation are plentiful [4J, even 
after the gauge freedom is quotiented out. The magnetic and scalar forces 
between static BPS monopoles cancel [5j , so there is great freedom in where 
the constituent monopoles can be located, and it has been rigorously estab- 
lished that there is a connected 4A^-dimensional moduli space of A^-monopole 
solutions. There is a deep mathematical theory, combining ideas of integrable 
systems, holomorphic bundles, Nahm data and spectral curves, and a mini- 
twistor space, giving insights into the Bogomolny equation and its monopole 
solutions [H |2]. The monopoles of charge N are in 1-1 correspondence with 
(based) rational maps of degree A^, from S"^ to S"^, known as Jarvis maps in 
this context [6]. By considering Jarvis maps, it becomes clear that no BPS 
monopole of charge greater than 1 can be exactly spherically symmetric, so 
a cluster of monopoles can be at most approximately spherically symmetric. 

The 1-monopole solution has a simple closed form, and much is known 
about 2- monopole solutions |T]. A 2-monopole solution has just one inter- 
esting parameter, the spatial separation of the pair. However, the the field 
values and energy density are complicated functions of spatial position and 
monopole separation. For charges higher than 2, only certain specially sym- 
metric monopoles are known precisely [2]. Even here, what is known are the 
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spectral curves and the rational maps associated with the monopole. The 
fields and energy distributions need to be calculated numerically. Examples 
are a few monopoles with platonic symmetry, with charges < 11. There 
is also a family of iV-monopole solutions, for all A^, which are regular poly- 
gons of unit charge monopoles [7j. When the radius of the regular iV-gon 
solution is minimised, the N monopoles merge, and the energy density is 
concentrated in a circular ring, or torus [SI El [ID] • The platonic monopoles, 
and the family of toroidal monopoles for all > 1, are among the most 
compact monopole solutions. 

If all A^ monopoles are very well separated from each other, then the so- 
lution is a superposition of A^ particle-like objects, each with little distortion. 
The 4N parameters are the A^ positions of the unit charge monopoles (3A^ 
parameters) and A^ phases which are required to define the gluing of their 
fields into space. However, there is little insight so far into the spatial struc- 
ture of a generic A^-monopole solution, where the monopoles form a relatively 
compact cluster. Monopoles cannot be very close together, as they get larger 
as they approach each other, so it is not clear how big such a cluster might 
be. 

The purpose of this paper is to shed some light on monopole clusters 
of large charge, and our method is to extend the analysis that was used 
by Bolognesi to discover and investigate a further class of BPS monopole 
solutions called magnetic bags [TT] . 

The existence of magnetic bags has not been rigorously proved, although 
it is likely they exist for all large A^. The bags are approximately spherical 
monopole solutions of charge A^, with the nonabelian structure [12] concen- 
trated on a thin spherical surface of radius approximately A^. Because of 
the surface structure, the fields are not exactly spherically symmetric. Bags 
can be distorted away from spherical, but we shall focus on those that are 
not. What is known precisely about magnetic bags are the abelianized gauge 
and Higgs fields that they produce. These resemble the electric field and the 
electrostatic potential arising from a charged conductor. 

The spherical magnetic bags are especially interesting, because, for each 
A^, they appear to be the most compact A^-monopole solutions that there 
are. We shall present new evidence for this in this paper. Our main pur- 
pose, however, is to generalise the construction of the magnetic bags to a 
much larger class of solutions that can be interpreted as spherical clusters 
of monopoles. Dilute clusters of A^ unit charge monopoles are very large, 
because the size of each constituent monopole depends on the Higgs field 
in the neighbourhood. If the background abelianized Higgs field is 0, then 
locally a monopole has radius 1/0. Our calculations below will give lower 
bounds on the size of monopole clusters based on this basic property of the 
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constituent monopoles. In a spherical cluster, is smallest near the centre, 
so here the monopoles are largest. 

In Section 2 we shall review the properties of the Bolognesi magnetic 
bag, and then in Section 3 present a multi-layer generalisation of the mag- 
netic bag, assuming spherical symmetry for simplicity. We shall establish an 
important constraint on the radial distribution of charge in such a structure. 
Various inequalities follow from this, and one of our conclusions will be that 
no multi-layer solution packs all monopoles into a smaller volume than 
the Bolognesi bag. In Section 4 we give a smoothed-out description of a 
multi-layer magnetic bag, which now becomes a cluster of monopoles with 
a quasi-smooth charge density. The density is a rather arbitrary function of 
distance from the centre, subject to one constraint. We discuss whether a 
Bolognesi type of bag persists in the inner part of the cluster. 

In Sections 5 and 6 we present examples of monopole clusters that have 
density profiles typical of two types of astrophysical object - the first is a 
planet, with a constant density inside and zero density outside; the second 
is a spherical galaxy, with a density falling with the inverse square of the 
distance from the centre [13]. We give these two types of cluster the names 
monopole planet and monopole galaxy. There are some analogies between 
monopole clusters subject to Yang-Mills-Higgs forces, and astrophysical ob- 
jects subject to gravitational forces, including black holes. We shall not 
simultaneously consider both types of force, although Bolognesi has recently 
considered gravitating magnetic bags [H]. 

In Section 7, which depends on more technical material about exact 
monopole solutions, we show, by calculating spectral radii, that the known 
monopoles with platonic symmetries are examples of approximately spherical 
magnetic bags, of small charge. We also show that the toroidal monopole of 
charge N is not as compact as the spherical magnetic bag of charge N. 

In Section 8 we present our conclusions, and summarise the analogy that 
appears to exist between monopole clusters and astrophysical objects, which 
extends the analogy between magnetic bags and black holes that has been 
noticed by Bolognesi |14] . 
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2 Magnetic Bags 



An interesting structure, discovered by Bolognesi and called a magnetic bag 
[TT] , is possible for a monopole with a large charge N. It is an abelian 
approximation to what is expected to be an exact solution of the SU{2) 
Bogomolny equation B = D$. The simplest version of the bag is spherically 
symmetric. The abelian fields are a scalar field 0, representing the magnitude 
of the Higgs field $, and a magnetic field b, the projection of the SU{2) 
magnetic field B onto the direction of $ in the SU{2) Lie algebra. They 
are related by b = V0, which is the abelianized Bogomolny equation. Away 
from the monopoles V ■ b = 0, so = 0; the monopoles are Coulomb- 
type sources for b. In the interior of the bag = 0, and in the exterior 
= 1 — N/r. Correspondingly, b vanishes in the interior and b = Nr/r^ 
in the exterior. At the bag surface, separating interior from exterior, is 
required to be continuous, so the bag is at radius R = N. Note that V0 and 
b are discontinuous across the bag, the discontinuity of the flux of b being 
47riV. The interpretation is that the bag carries magnetic charge N, spread 
uniformly over the surface, and that the interior of the bag is a magnetic 
conductor, carrying no charge. 

It is very plausible that there is an exact solution, approximated by the 
bag, which has abelianized fields as described above except in a thin spher- 
ical wall at radius A^, where the fields are fully nonabelian. The bag wall 
is of thickness approximately \/N, and is a slightly curved version of the 
fiat monopole wall solution found by Ward [15]. That solution has a two- 
dimensional lattice periodicity, which is why the bag cannot be exactly spher- 
ically symmetric. The spherical bag must not only curve the wall, but also 
introduce defects in the two-dimensional lattice structure. Nevertheless, the 
lattice structure probably survives locally. Lee and Weinberg have described 
in considerable detail some of the topological features of the nonabelian fields 
of the bag [12]. Nahm data associated with magnetic bags of large charge 
have been discussed by Harland [16J. 

Some evidence for the existence of bags comes from the platonic monopoles 
[3 [TTl [18] . These are exact solutions of the Bogomolny equation of low charge 
which resemble bags, in that they are hollow platonic polyhedra, with energy 
and nonabelian field structures concentrated on their surfaces. So far, their 
radii have not been compared with the Bolognesi estimate R = N. Using 
the spectral curves of these monopoles, we calculate their radii in Section 7 
here, and show that they are close to the Bolognesi value. As a definition of 
size we use the spectral radius [1]. 

Perhaps the most surprising property of a magnetic bag is its large size. 
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Since the radius is A^, its area is a multiple of A^^ and its volume a multiple 
of A^^. One might expect BPS monopoles to be a model of rather normal 
matter, where particles fit into a volume of order A^, but this is not the 
case. The magnetic bag is much larger. Black holes have a similarly exotic 
size scaling. The radius of the horizon of a Schwarzschild black hole of mass 
M is 2M in natural units, and the black hole area is proportional to M^. 
A magnetic bag is therefore, to some extent, a gauge theory analogue of a 
black hole. 

It is worth recalling a few more properties of the magnetic bag, and some 
variants of the bag, pointed out by Bolognesi, and by Lee and Weinberg. 
First, the bag can be deformed from spherical, because of the large number 
of dimensions of the moduli space of A^ monopoles. However, we will only 
consider spherical distributions of monopoles here. Next, there is a rather 
trivial rescaling of the Bolognesi bag if the vacuum expectation value of 
the Higgs field $ is rather than the special value v = 1. One now has 
= v — N/r outside the bag, and = inside, so the bag radius is i? = N/v. 
Its thickness is a/ZV/u [T2]. The argument for this is that the bag area is 
N'^/v'^ (dropping constants of order 1), so the area per unit charge monopole 
is N/v'^. The characteristic length associated with a monopole in the bag 
surface (which can be thought of as the monopole separation within the 
surface) is therefore the square root of this, a/ZV /v, and since it is the only 
length scale in the bag, the bag thickness must be the same. We will, in fact, 
always suppose that v = 1 and impose the boundary condition that — )■ 1 
as r — )■ oo. However, the bag solution with other values of v will be locally 
relevant. Note that for any v the bag thickness is R/\/N, where R is the 
radius and A^ the charge. 

Lastly, we recall that there is a variant bag of charge A^ where does not 
vanish inside, but instead = c inside, with < c < 1. Outside, = 1 — N/r 
as usual, so the bag radius is R = N/ (1 — c), making continuous across the 
bag surface. This bag radius is larger than the radius of the Bolognesi bag, 
and if c is not very small, then the bag breaks up into constituent, unit charge 
monopoles. To estimate the critical value of c where break-up occurs, note 
that since the Higgs field is = c in the bag vicinity, unit charge monopoles 
have a radius 1/c. The bag has broken up into constituents if the monopole 
separation within the bag, Vn, is bigger than this. The critical value of c is 
therefore 1/y/N, which is much less than 1 if A^ is large. The critical value of 
c is reached when the bag has expanded to a radius approximately A^ + 
If the bag is larger than this it should be regarded as a spherical layer of 
separated, unit charge monopoles, rather than as a Bolognesi bag. 
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3 Mult i- layer Magnetic Bags 



Using the abelianized Bogomolny equation b — V0, it is straightforward to 
calculate the Higgs field and magnetic field b of a multi-layer structure of 
several magnetic bags with spherical symmetry. Let us suppose that there 
are bags with positive integer charges {Ni, N2, . . . , Nk} on spheres of radii 
{Ri, R2, • • • , Rk}, where < i?i < i?2 < • • • < Rk- The total charge is 

N — Ni + N2-\ h Nk- Let b denote the magnitude of the radial magnetic 

field b. b{r) has discontinuities at each bag. Between R^ and Rk+i, 

Ni + N2 + --- + Nk 

= -2 ' (1) 

the magnetic analogue of the Coulomb field due to the enclosed charge. Also 
b{r) = for r < i?i and b{r) = N/r'^ for r > Rk- The Higgs field 0(r) 
satisfies 

^ = (2) 

If we integrate out from r = 0, and assume 0(0) = 0, then we find that 
4>{r) = for r < Ri, 

^M.^ + ^ + ... + ^_^XL±*+:i±JYi (3) 

Ri R2 Rk T 

between Rk and Rk+ii and 

' R^^ R2 Rk r ^ ' 

for r > Rk- (p is continuous at each bag radius, but d(f)/dr has a discontinuity. 
The boundary condition — > 1 as r — > 00 requires that 

ril 1x2 -tCk 

This is the key constraint on multi-layer bags and their distribution of charges. 
This constraint is the generalisation of the result that for a single Bolognesi 
bag, R — N. Various inequalities follow. Among them we have 

-^<h (6) 

so each bag is larger than it would be if it had the same charge but were 
alone. More interesting is to consider the total charge N. Replacing each 



6 



radius Rk by Ri we increase the sum on the left hand side of eq.dS]) to a 
value greater than 1. But the sum is now N/Ri, so Ri < N. Similarly, 
replacing Rk by Rk, we find Rk > N. Thus, compared with a Bolognesi bag 
of charge N, the multi-layer structure has a smaller innermost bag, and a 
larger outermost bag. The structure as a whole is larger than the Bolognesi 
bag. Indeed, this is strong evidence that the smallest radius into which one 
can fit monopoles is the Bolognesi radius R = N. 

A very similar solution with </)(r) = c for r < Ri (with < c < 1) exists 
when 

^ + ^ + --- + ^ = l-c- (7) 

Iti 112 

Here Rk > N/ (1 — c), so the solution is larger than before. 

Let us now consider the type of structure each bag layer has. For the rest 
of this section we assume c = 0. The innermost layer is a bag of Bolognesi 
type. It has charge Ni and radius Ri, the Higgs field is = inside and 
0(r) = Ni/ Rl — Ni/r outside (as far as the second layer). This is just as 
expected for a Bolognesi bag of charge A^i and effective vacuum expectation 
value V = Ni/Ri. The thickness of this bag is \/Wi/v = Ri/\/7^, which is 
small compared to Ri provided A^^i is large. The second layer needs to be 
further from the first than this thickness, for the multi-layer point of view to 
be justified. 

Let us therefore assume that the separation R2 — Ri is significantly larger 
than Ri/^/Ni, but significantly less than Ri. Let us also assume that the 
charges in the first and second layer are similar. Within the second layer 
the monopole separation is R2/y/7^, which is comparable to the separation 
within the first layer. Just inside the second layer the Higgs field is 

^ Ni iVi NiiR2-Ri) 



Rl i?2 R1R2 

which, by our assumptions, is considerably larger than both y/N\/R2 and 
y/N2/R2- We estimate the monopole size in the second layer as 1/0 = 
R1R2/ {Ni[R2 — Rl)), and this is considerably less than the monopole sepa- 
ration within the layer, i?2/\/^- The second layer is therefore not a Bolog- 
nesi bag, but instead a spherical arrangement of well separated unit charge 
monopoles. Since increases as one moves out to subsequent layers, the char- 
acteristic size of a monopole decreases. On the other hand, the monopole 
separation within each layer increases as one moves out, provided the charges 
in all the layers are similar. It follows that all the outer layers consist of well 
separated monopoles, rather than a bag of Bolognesi type. The entire so- 
lution consists of an inner layer which is a Bolognesi bag, surrounded by 
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a layered cluster of well separated monopoles. These monopoles are not of 
equal size. They are largest near the centre, where the Higgs field is small, 
and decrease in size further out, with monopoles in the outer layers having 
approximately unit size. 

A multi-layer bag solution may have a Bolognesi type of bag as one of 
the intermediate layers. But this only happens if the distribution of charges 
is very uneven, with the constraint satisfied by having one ratio Nk/R^ 
very close to 1, and all the others much less than 1. The kth layer is then 
a Bolognesi bag. It is still possible for one of the layers closer to the centre 
to be a Bolognesi bag too, but its charge must be very much smaller. We 
can verify this by considering a structure with just two layers, each with a 
large charge (so that the square root of the charge is small compared with 
the charge itself). If we fix their radii at Ri and R2 (with ratio of order 1), 
and satisfy the constraint 

by setting Ni/Ri = e and N2/R2 = 1 — e with e small, then the inner layer 
is a Bolognesi bag automatically, and the outer layer is one too, provided 
is sufficiently small at r = i?2- This occurs if e < I/a/A^ . Since N2 is 
approximately R2, and Ni = eRi < eR2 — BN2, Ni needs to be less than 
v^iV2. As anticipated, the inner charge is much smaller than the outer charge. 

Our main conclusion is that multi-layer bag structures are possible, but 
the layers are not generally of the Bolognesi bag type. Typically, the in- 
nermost layer is a Bolognesi bag, and it is surrounded by spherical layers 
of isolated unit charge monopoles. We shall next consider the continuum 
limit of these structures. There will be a spherically symmetric density of 
magnetic charge located between an inner and outer radius. At the inner 
radius we shall find there is often a Bolognesi bag, and outside this, a gas of 
isolated monopoles. 



4 A Smooth, Spherical Cluster of Monopoles 

For a spherical cluster of BPS monopoles, of very large charge, the continuum 
limit arises when the layers discussed above merge into a charge density that 
is a quasi-continuous function of the radius. It is convenient to work with the 
radial density q{r), defined so that the charge in the shell between radius r 
and radius r+dr is q{r)dr. For monopoles satisfying the Bogomolny equation. 
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g(r) > everywhere. The charge Q{R) interior to radius R is 

Q{R)= rq{r)dr, (10) 



^0 

and the magnetic field strength is b{R) = Q{R)/R^. We shall assume that 
the total charge, Q = limji^aoQiR) = N, is finite. 

The constraint becomes, in the continuum limit, 

'-^dr = l. (11) 
r 

This is in the case that the Higgs field runs from at the origin to 1 at 
infinity. More generally, if (p has a positive value c at the origin, 

^c/r = l-c. (12) 
r 

These integrals automatically converge at infinity if the total charge is finite, 
but will fail to converge at the origin unless the density q tends to zero 
sufficiently fast there. Having a finite, non-zero volume density of charge 
at the origin makes g(r) proportional to for small r, and this is one way 
to ensure convergence. We shall also consider q strictly vanishing inside an 
interior radius Ri^. It also makes sense to assume that q vanishes outside 
some outer radius -Rout- This is because monopoles each have at least one unit 
of charge, so that a smoothly decaying charge density stretching to infinity 
is not possible. 

With the outer radius assumed, the total charge is 

/"Rout 

Q= / q{r)dr, (13) 



and the constraint ffT^ is 







r 



dr = l-c. (14) 



Replacing r by i?out in the denominator gives the inequality Q /Rout < 1 — c 
(assuming q not entirely concentrated at -Rout), so the ball that contains all 
the charge has radius larger than Q/(l — c), the radius of the one-layer bag 
for the same total charge. 

We will need the formula for the Higgs field 0. By analogy with the 
multi-layer case, or by direct integration, we find 

riw*^5M+,. (15) 

.In f R 
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This is analogous to the Newtonian gravitational potential for a spherical 
mass distribution, shifted so that = 1 at infinity. Using the constraint 
(fT¥|) . we have the alternative expression 



cf>{R) = 1 



/ 



q{r) 



dr — 



Q{R) 



(16) 



r 



R 



(f){R) obeys the spherically symmetric version of the Poisson equation 



and b{R) = d(f)/dR as before. 

For generic densities g(r), we anticipate that if c = a Bolognesi bag 
could occur at some inner radius, with a gas of isolated monopoles outside. 
For sufficiently non- uniform q{r) there could be Bolognesi bags at different 
radii, as we saw with the multi-layer structures, but we will not seek these out 
again. In the next sections, we consider two specific densities, corresponding 
to what we think of as a monopole planet and a monopole galaxy. The 
monopole galaxy has a Bolognesi bag at its inner radius. 



A (small) planet is a spherical body where gravity plays a limited role. If the 
material is uniform, the density is uniform, dropping suddenly to zero at the 
planet's outer radius. By a monopole planet, we mean a spherical cluster of 
monopoles with a radial density g(r) = 47rpr^ between r = and r = -Rout, 
with p, the charge density per unit volume, a positive constant. The density 
vanishes outside. As the volume per unit charge is we estimate the 
nearest neighbour monopole separation as p^^/^. Here, and in the rest of 
this and the next section, factors of order 1 are often ignored. 
The total charge of the planet is 




(17) 



5 Monopole Planets 



Q = -jpRl 



(18) 



and the constraint f|T^ simplifies here to 



27rpi?2 



1 -c. 



(19) 
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Eliminating p in favour of Q we find 



3 



Q. 



(20) 



'Out 



2(1 -c) 



This is always greater than the Bolognesi radius R = Q. 

One might imagine that for given charge Q, -Rout is smallest when c = 0, 
but this choice of c is not consistent. Close to the origin (p would be so 
small that the monopole sizes would exceed their separation. To restore 
consistency we would need to dilute the monopole density near the origin, 
creating a hole. We do not want to do this here, but will do so in the 
next section. Instead, let us keep the density uniform, but find the minimal 
consistent value of c. This is the value c = p^^^, for which the monopole 
size close to the origin, which is 1/c, matches the monopole separation p~^/^. 
Again we can eliminate p in favour of Q. We find that c = 0(Q-2/3) if 
Q is large. So the monopoles close to the centre of the planet have size of 
order Q'^/^ (and hence volume Q^), much larger than the monopoles near the 
outside, which have unit size. From f l2U]) it follows that this densest possible 
monopole planet has size 

near the centre but not zero, so there is no magnetic bag in the centre. A 
less dense monopole planet of the same charge arises if c has a larger value. 



6 Monopole Galaxies 

The remarkable feature of many spherical and almost spherical galaxies is 
that their density is highly non-uniform, decaying as the inverse square of 
distance from the centre ^I3j . This density profile is called a singular isother- 
mal sphere. We shall investigate here the structure of a BPS monopole of 
large charge with such a density profile. If the charge density per unit volume 
decays as then the radial density g(r) is a constant, go- Let us suppose 
that g(r) = go in the radial range i?in < r < i?out and vanishes outside this 
range. The total charge Q is go (-Rout — -Rin)- We want a finite charge, so -Rout 
is finite. If -Rin ^ -Rout then the total charge is approximately go-Rout- 
The usual constraint on the radial density becomes 




(21) 



that is. 



go (log/?, 



'Out 



log -Rin) = 1 - C . 



(22) 



11 



This constraint forces i?in to have a positive value. Having an inner cut-off 
at a positive radius is the same as in galaxy modelling, and avoids a singular 
density at the origin. The interesting case is where c = 0. The Higgs field 
then vanishes for r < Ri^, and we anticipate that there can be a Bolognesi 
bag in the neighbourhood of R^, surrounded by a more dilute gas of well 
separated monopoles. Note that go is not very different from 1 for a large 
range of ratios Rout/Rin- 

We need the Higgs field (p between _Rin and -Rout- The expression fflGjl 
reduces to 

(/.(i?) = 1 - go (log i?out - log i?) - go f 1 - %1 • (23) 



R ^ 

For R near i?out this has the characteristic logarithmic dependence on R 
familiar from the Newtonian gravitational potential in galaxy models. Near 
i?in there is a quadratic dependence on R — Ri^, and the slope of (to leading 
order) is 

^^f (24) 

Now recall that for the multi- layer structures discussed in section 3, with 
= in the interior, the innermost layer is usually a Bolognesi bag. The 
charge A^i and radius Ri are arbitrary, but they are related to the slope 
of just outside the bag, which is Ni/Rf. The continuum distribution of 
charge density also allows a Bolognesi bag by this criterion, as the following 
argument shows. Suppose the bag has thickness /, small compared with 
so the bag lies between Ri^ and Rin + I- The charge in the bag is q^l and 
the bag radius is -Rin. The expected slope of the Higgs field just outside the 
bag is qol/Rf^. But this agrees with the slope fl2^ . So a Bolognesi bag is 
present, but this calculation doesn't determine its thickness, nor its charge. 
However there is a further consideration, as we have specified the charge 
density throughout the monopole galaxy, and in particular, near the inner 
radius. The volume per unit charge here is Rf^/qo (dropping a in factor), so 
the separation of monopoles is 



p2/3 
R-m 

1/3 
% 



(25) 



If we identify s with -Rin/ y/qol, the separation of monopoles within a Bolognesi 
bag of radius -Rjn and charge go/, we deduce that I = Rfl^ /ql^^ = s. The 
thickness of the bag I is the same as the monopole separation within the bag, 
which is what we expect. The charge of the bag is Qbag = {QoRinY^^- 
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So we have now consistently determined the parameters of the Bolognesi 
bag at the centre of a monopole galaxy. In terms of the total charge Q and 
the inner and outer radii of the monopole galaxy, the bag charge is 



2/3 



(26) 



using the approximation that Q = go-Rout- Qbag is very much smaller than 
Q, so almost all monopoles form a dilute gas outside the bag. Similarly, the 
thickness is 



Therefore, a monopole galaxy has a feature rather analogous to the black 
hole at the centre of a real galaxy, namely an inner magnetic bag where 
monopoles have coalesced, carrying a small fraction of the total charge. Given 
the bag radius, and the outer radius and total charge of the monopole galaxy, 
the charge of the bag is determined. This contrasts with the black hole in a 
galaxy, whose mass is not theoretically determined. 

We can have some fun with the monopole arithmetic. Suppose a monopole 
galaxy has charge 10^^, and the ratio of outer to inner radius is 10"^. Then 
the magnetic bag at the centre has charge 10^. These numbers resemble the 
analogous numbers for a real galaxy and its black hole (with charge identified 
with number of solar masses). So the monopole equations, which forbid (j) to 
be negative, provide a useful model for regularising the density singularity 
inside a galaxy, as well as providing an interesting analogy for understanding 
why central black holes occur in galaxies. 



7 Spectral Radii of Monopoles 

This section depends on the theory of monopole spectral curves, and the 
details, other than the next paragraph, may be skipped. 

There are BPS monopole solutions of charges = 3, 4, 5, 7 with the sym- 
metries and shapes of the platonic solids. They provide small scale models 
of the Bolognesi bag that is expected at large N. Here we use the known 
spectral curve of each of these monopoles to calculate its spectral radius D^, 
which is a precisely defined estimate of the monopole size, and compare with 
the Bolognesi estimate of the radius, Rn- We will show them to be similar. 
We include the = 1 and toroidal N — 2 monopoles in the calculations. 




(27) 
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and also make some remarks about the sequence of axisymmetric, toroidal 
monopoles that exist for all A^. 

Any BPS monopole is characterised by its spectral curve. This is ob- 
tained easily from the Nahm data of the monopole, or equivalently, using 
the normalisable solutions of the Hitchin equation in the monopole back- 
ground. For details see, e.g. [U The spectral curve is a complex curve 
lying in the space of all straight lines in R'^ (which is the mini-twistor space 
TP^, a non-compact complex surface). The complex coordinates of TP^ are 
parametrising the direction of a straight line (obtained by usual stereo- 
graphic projection), and rj which parametrises the intersection point of the 
line with a (complex) plane orthogonal to the direction ^. Given a line with 
coordinates ^ and r], the points {xi,X2,X3) G R'^ that lie on it are those 
satisfying 

r/ - (x2 - ixi) + 2x3^ + {X2 + ixi)^'^ = . (28) 

Important for us is to note that the distance between the line and the origin 
in R^ (the distance to the closest point on the line) is 

If ?7 = then the line passes through the origin. 

The spectral curve of a given monopole of charge is an algebraic curve 
in TP^ of the form 

PNiv,O = 0, (30) 

where Pn is a polynomial of degree N inr] and (generically) of degree 2N in ^. 
For fixed direction ^, there are solutions rj. So the monopole has straight 
lines associated with it in this direction (and in every other direction). These 
are called spectral lines. Heuristically, spectral lines pass through regions 
where the energy density of the monopole is concentrated and the fields are 
fully non-abelian. In particular, for A^ well separated unit charge monopoles, 
these lines pass through all the monopole cores. 

The spectral curve is compact. This implies that the distance d defined 
above has a finite maximum value D as and rj vary over the spectral curve. 
Let us call D the spectral radius of the monopole. There is a notion of a 
centred monopole (with centre of mass at the origin). The spectral radius 
is most useful when applied to a centred monopole, and is invariant under 
rotations of the monopole. 

Let us now evaluate the spectral radii of the most symmetric, centred 
monopoles of low charge. We consider the basic A^ = 1, 0(3)-symmetric 
monopole, the toroidal N = 2 monopole, and the bag-like monopoles of 
charges A^ = 3,4,5,7 with tetrahedral, cubic, octahedral and dodecahedral 
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form [3 [T71 HB]- There is also a monopole of charge = 11 with icosahedral 
form [12]. They are unique for these charges, and no monopoles of smaller 
charge have these platonic symmetries. 

The spectral curves of all these have been determined, except the = 11 
example. They are, respect iveljH, 

V = 0, (31) 
v' + ^e = 0, (32) 

r (if r (if 

+ ^^l^^ae - I) = 0, (33) 
48v37r3/2 

2567r2 



v' + ^^ie + Ue + l)v = 0, (35) 



v'+ 'q^J' ae'+ne-i)v = o. (36) 

Use of the identity r(l/6) = {3/7iY/'^2-'^/'^iT{l/3)Y gives alternative expres- 
sions for the constants above [21]. In the cases N > 2, the polynomial /C(0 
that appears here is a Klein polynomial associated with the relevant platonic 
solid [22]. This is a polynomial that vanishes in the directions of all the face 
centres, or all the edge centres, or all the vertices of the solid. Note that we 
are following the conventions of [2J in this section, where the asymptotic be- 
haviour of the Higgs field is = 1 — N/2r. The Bolognesi radius is therefore 
Rn = N/2. 

For given ^ it is easy to find the roots r) for each of these equations. One 
of the roots is r/ = in cases fl35p and fl36p . The non-zero roots all have 
the same magnitude [t^I^. Evaluating = \ri\^/{l + |^p) therefore gives the 
maximum distance from the origin of spectral lines in the direction ^. is 
a function with the symmetry of the monopole. The spectral radius, which 
we denote hj Dn for these special monopoles, is the maximum value of d^ 
over all values of ^. 

For N = 1, d^ = for all ^, so = 0. For = 2 



^Included in equations (p4l) and ((35)) is a factor 16 correction to equations in [71 [T8]. 
see [201. 
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and this is maximal on the equator |^| = 1. So 

= ^ ~ 0.785 . (38) 

For higher N it is not always easy to find hy elementary calculus. A 
reasonable conjecture is that it occurs in a direction with a simple meaning 
for the relevant platonic solid. Let V, E and F denote the sets of directions ^ 
corresponding to the vertices, edges and faces of the solid. Let ny, ue and np 
denote the number of these elements. By Euler's theorem, riy — n-E + n-p = 2. 
Because |?7|^ has no branch crossing as ^ varies, is a smooth Morse function 
on the 2-sphere. Denote the number of its local maxima, saddle points and 
minima by nmax, risaddie and n^in- By the Morse relations, ^max - rigaddie + 
n-min — 2. For a simple function with the symmetry of the platonic solid, 
two possible arrangements of the stationary points are therefore that the 
maxima occur at V(or F), the saddle points at E and the minima at F(or 
V) . Exchanging V and F means replacing the sohd by its dual. 

We will assume that for the spectral curves above, is such a simple 
Morse function. In each case the minimum value of d^ is zero, the minima 
occuring at the roots of the Klein polynomial /C(0 appearing in the spectral 
curve equation. For = 4, A^(^) = + 14^"^ + 1 is the vertex polynomial 
of the cube, so the minima are at the vertices, and the maxima at the face 
centres of the cube. For N — 5, /C(^) is the same, so the conclusion is the 
same (although we can now say the minima arc on the faces of an octahedron 
and the maxima at the vertices). For N = 7, /C(0 = ^{^^^ + 11^^ — 1) is 
the face polynomial of the dodecahedron (with one root at infinity), so the 
maxima of d^ are at the vertices of the dodecahedron. 

The case = 3 is more tricky. While the monopole itself has tetrahedral 
symmetry, and is not invariant under an inversion, the function d^ has this 
inversion symmetry, because if a line in IR^ belongs to the spectral curve, 
then so does the line obtained by reversing its direction. Inversion sends ^ to 
its antipode — The distances d for this pair are the same. The function 
d^ therefore has not just tetrahedral symmetry, but cubic symmetry. For 
N — 3, JC{^) — ^({^ — 1) is the edge polynomial of the tetrahedron, so the 
edge centres are where d^ has its minima. However, these points are also the 
face centres of the cube whose vertices combine those of the tetrahedron and 
its dual. The maxima of d^ occur on the vertices of this cube, that is, on 
both the faces and vertices of the initial tetrahedron. (The 12 saddle points 
are at the edge centres of the cube, and are distributed three per face of the 
tetrahedron.) Having identified the locations of the maxima of d^ for each of 
the platonic monopoles we can now calculate their spectral radii. 
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For N — 3, one vertex of the cube just described is at 

e=^(V3-l), (39) 

this being a root of the Klein polynomial + 14^"^ + 1. Solving for |?7|^ from 
the N — 3 spectral curve gives 

^ 24/3(37r)V2 l+|^|2 ' ^^^^ 
and evaluating this with the above value of ^, we find 

r (i) r (i) 2^/6 

' 2V3(37r)V2 3V2 - ^^^^ 

It is easy to verify directly that this is the maximum value of d^. Firstly, for 
given 1^1, ^'^ should be negative, which determines the possible arguments of 
^, and then the value of |,^| which maximises the remaining real expression 
can be found by differentiation. 
For AT = 4 

= 1^1024;^ ) TTp ■ ^^^^ 

The maxima of d^ are at the face centres of the cube, the roots of ^(^^ — 1), 
one of which is at ^ = 0. Hence 



1/4 

D4 = I 7noV^ I - l-'^25. (43) 



\^10247r2 j 

One can verify that d^ has the same value at other face centres, e.g. at ^ = 1, 
and it has a lower value at the edge centres, e.g. at ^ = -\/2 — 1. Again, D^^ 
can be directly verified to be the maximum value of d^. The maximum occurs 
for ^'^ real and non-negative, so we can choose ^ real. Then we can find the 
stationary points of {^^ + 14^"^ + 1)^/"^/ (1 + ^^) by differentiation, and identify 
the maxima at ^ = 0, ±1 (and ^ = 00). 

For N — 5 the geometry is the same, but the spectral curve equation has 
an extra factor of 4 (and an additional overall factor of 77). This means that 
the non-zero roots 77 are multiplied by ^/2, so D5 = \f2D/^. Therefore 

D5 ~ X 1.725 ~ 2.440 . (44) 

For N = 7 we need a vertex of the dodecahedron, a root of the Klein 
polynomial _ 228^^^ + 494^^° + 228^^ + 1. We have found numerically 
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that one of the roots is ^ = -0.3382612 with = -0.00442854 [23j. Solving 
the spectral curve equation gives 

, r(l)r(|) m'' + lie - 

27rV2 l+|^|2 ' 

and evaluating this for the numerically obtained ^ gives 

L>7^ 3.176. (46) 

Here we have not shown directly that this is the maximum value of d^. 

Our combined results are that for N = 1,2,3,4,5,7, the values of the 
spectral radius are 0,0.785,1.249,1.725,2.440,3.176. For comparison 
the Bolognesi radius is Rn = N/2, with values 0.5, 1.0, 1.5, 2.0, 2.5, 3.5. The 
ratios D^/Rn are 0,0.785,0.833,0.863,0.976,0.907, and these appear to be 
approaching 1 as increases. These calculations therefore confirm that 
the platonic monopoles are models for the Bolognesi magnetic bag at small 
charge. 

Finally, we consider the spectral radii of the sequence of toroidal monopoles 
that exist for all A^. Using the spectral curves obtained by Hitchin [10], one 
sees that the largest roots rj for any C, are 

- 1 

^ = ±^— vre, (47) 

so |?7|g = (A^ — l)7r|,^|/2. It follows that the maximum value of = 1^71^/(1 + 
|,^P) is (A^ — l)7r/4, and this is the spectral radius of the charge A^ toroidal 
monopole. Bolognesi has estimated this using an abelian magnetic disc ap- 
proximation, finding a radius Nir/A (rescaling by 2, as elsewhere in this sec- 
tion) [11] . These radii agree asymptotically, the spectral radius being slightly 
smaller as for the platonic monopoles. 

We see that for any A^, the toroidal monopole in the magnetic disc ap- 
proximation is larger than the spherical magnetic bag by a factor 7r/2, which 
is further evidence that the spherical magnetic bag is the most compact 
monopole. We can also verify that for A^ = 3,4,5,7 the spectral radius of 
the toroidal monopole is larger than that of the platonic monopole. 



8 Conclusions 



Magnetic bags are clearly an important type of BPS monopole solution, even 
if their existence has not been rigorously established. The approximately 
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spherical magnetic bag, for each large charge N, is probably the most com- 
pact monopole within the whole 4 A^- dimensional moduli space of solutions, 
and we have found new evidence for this. 

A large class of new monopole solutions has been found, using the same 
abelian approximation that gave the magnetic bag solution. These are spher- 
ical clusters of monopoles, and they are all larger than the Bolognesi bag of 
the same charge. Some of these clusters have bag structures in their interior, 
and it would be interesting to investigate in more detail how the bag joins 
on to the gas of more isolated monopoles outside the bag. 

We have shown here that the known platonic monopoles with = 
3,4,5,7 form a sequence that look bag-like, and that with increasing A^ 
their spectral radii approach the asymptotic value of the radius predicted 
for a magnetic bag, R = N. There should be further, exact, bag-like so- 
lutions with platonic symmetry. The first of these has A^ = 11. There is 
then a sequence of rational maps, with icosahedral symmetry and various 
degrees from 17 up to 97 and beyond, which have been used to construct 
bag-like Skyrmions . These could produce bag- like monopoles too, if one 
interpreted these rational maps as the Jarvis maps of monopoles, and imple- 
mented the monopole construction of loannidou and Sutcliffe ^25]. Less clear 
is the characterisation of the Jarvis maps which correspond to the spherical 
monopole clusters of variable density that we have discussed. A spherical 
cluster would probably arise from a Jarvis map whose zeros and poles were 
distributed randomly but evenly over 5*^. However, it is difficult to construct 
a monopole of large charge from a rational map of large degree, as it is still 
necessary to solve a partial differential equation numerically. 

We have interpreted the monopole cluster solutions found here in astro- 
physical terms. This was motivated by the apparently good analogy between 
a magnetic bag of charge A^ and a black hole of mass M [H] . The bag radius 
is R = N and (in units where Newton's constant and the speed of light are 
1) the horizon radius of the black hole is R = 2M. 

For the monopole clusters, we have seen that (p runs either from at the 
centre to 1 at infinity, or from c at the centre to 1 at infinity, with c > 0. 
Moreover, if = at the centre, then there is a Bolognesi bag there. These 
possibilities are analogous to the result that the Newtonian gravitational 
potential of an astrophysical cluster (in our units) runs from some negative 
value no less than —1/2 at the centre to at infinity, and if the central value 
is —1/2, then light cannot escape, and there is a central black hole. We have 
seen that the objects we call monopole planets have c > at the centre, 
so these are like real planets or stars that have no black hole inside. The 
monopole galaxies, on the other hand, can have c = 0. The inner part of 
such a monopole galaxy is a magnetic bag, and this is analogous to a black 
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hole at the centre of a real galaxy. 

The conclusion seems to be that clusters of BPS monopoles of large mag- 
netic charge provide a new, non-relativistic analogy for both black holes and 
less dense astrophysical objects. The analogy does, of course, have limita- 
tions. In particular, the monopole solutions considered here are all static, 
whereas clusters of stars could only be approximately stationary if the indi- 
vidual objects in the cluster were orbiting the centre. Also, unfortunately, 
monopoles are not physically observable, although they can be studied math- 
ematically, and numerical solutions presented visually [2]. 

In summary, the analogy rests on three points. First, both BPS monopoles 
and astrophysical objects can be approximately described by an abelianized 
scalar potential satisfying a Poisson equation. In the monopole case (f) is 
the abelianized Higgs field; in the astrophysical case it is the Newtonian grav- 
itational potential in natural units. In both cases, has to be course-grained, 
to smooth out the sharp dips due to unit charge monopoles in monopole clus- 
ters, or stars within galaxies. Second, the potential cannot be less than a 
critical value: in the monopole clusters, and —1/2 in astrophysical objects. 
In the monopole case this is not so much because the true Higgs field $ has 
minimum magnitude zero at the centre of monopoles, but rather because 
determines the local size of a unit charge monopole, which has to remain 
finite; in the gravitational case it is because the Newtonian approximation 
breaks down where reaches —1/2, as light cannot escape from such a region, 
and the full machinery of Einstein's equations is needed. Third, the region 
where has its critical value is a region of very compact matter, and ex- 
hibits an exotic size scaling, with its radius rather than volume proportional 
to the amount of matter. For monopoles, such a region is the interior of a 
Bolognesi magnetic bag, the most compact arrangement of magnetic charge. 
The radius of the bag is equal to the bag's magnetic charge. For gravitating 
matter, the region is the horizon of a black hole (with the Newtonian picture 
invalid inside), whose radius is twice the black hole's mass. 
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